home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_her2.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  5.0 KB  |  196 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_her2 (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int order = 101;
  13.    int uplo = 121;
  14.    int N = 1;
  15.    int lda = 1;
  16.    float alpha[2] = {-1.0f, 0.0f};
  17.    float A[] = { -0.821f, 0.954f };
  18.    float X[] = { 0.532f, 0.802f };
  19.    int incX = -1;
  20.    float Y[] = { 0.016f, -0.334f };
  21.    int incY = -1;
  22.    float A_expected[] = { -0.302288f, 0.0f };
  23.    cblas_cher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
  24.    {
  25.      int i;
  26.      for (i = 0; i < 1; i++) {
  27.        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher2(case 1450) real");
  28.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher2(case 1450) imag");
  29.      };
  30.    };
  31.   };
  32.  
  33.  
  34.   {
  35.    int order = 101;
  36.    int uplo = 122;
  37.    int N = 1;
  38.    int lda = 1;
  39.    float alpha[2] = {-1.0f, 0.0f};
  40.    float A[] = { -0.821f, 0.954f };
  41.    float X[] = { 0.532f, 0.802f };
  42.    int incX = -1;
  43.    float Y[] = { 0.016f, -0.334f };
  44.    int incY = -1;
  45.    float A_expected[] = { -0.302288f, 0.0f };
  46.    cblas_cher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
  47.    {
  48.      int i;
  49.      for (i = 0; i < 1; i++) {
  50.        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher2(case 1451) real");
  51.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher2(case 1451) imag");
  52.      };
  53.    };
  54.   };
  55.  
  56.  
  57.   {
  58.    int order = 102;
  59.    int uplo = 121;
  60.    int N = 1;
  61.    int lda = 1;
  62.    float alpha[2] = {-1.0f, 0.0f};
  63.    float A[] = { -0.821f, 0.954f };
  64.    float X[] = { 0.532f, 0.802f };
  65.    int incX = -1;
  66.    float Y[] = { 0.016f, -0.334f };
  67.    int incY = -1;
  68.    float A_expected[] = { -0.302288f, 0.0f };
  69.    cblas_cher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
  70.    {
  71.      int i;
  72.      for (i = 0; i < 1; i++) {
  73.        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher2(case 1452) real");
  74.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher2(case 1452) imag");
  75.      };
  76.    };
  77.   };
  78.  
  79.  
  80.   {
  81.    int order = 102;
  82.    int uplo = 122;
  83.    int N = 1;
  84.    int lda = 1;
  85.    float alpha[2] = {-1.0f, 0.0f};
  86.    float A[] = { -0.821f, 0.954f };
  87.    float X[] = { 0.532f, 0.802f };
  88.    int incX = -1;
  89.    float Y[] = { 0.016f, -0.334f };
  90.    int incY = -1;
  91.    float A_expected[] = { -0.302288f, 0.0f };
  92.    cblas_cher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
  93.    {
  94.      int i;
  95.      for (i = 0; i < 1; i++) {
  96.        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher2(case 1453) real");
  97.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher2(case 1453) imag");
  98.      };
  99.    };
  100.   };
  101.  
  102.  
  103.   {
  104.    int order = 101;
  105.    int uplo = 121;
  106.    int N = 1;
  107.    int lda = 1;
  108.    double alpha[2] = {-0.3, 0.1};
  109.    double A[] = { -0.334, 0.286 };
  110.    double X[] = { -0.14, -0.135 };
  111.    int incX = -1;
  112.    double Y[] = { 0.455, 0.358 };
  113.    int incY = -1;
  114.    double A_expected[] = { -0.264521, 0.0 };
  115.    cblas_zher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
  116.    {
  117.      int i;
  118.      for (i = 0; i < 1; i++) {
  119.        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher2(case 1454) real");
  120.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher2(case 1454) imag");
  121.      };
  122.    };
  123.   };
  124.  
  125.  
  126.   {
  127.    int order = 101;
  128.    int uplo = 122;
  129.    int N = 1;
  130.    int lda = 1;
  131.    double alpha[2] = {-0.3, 0.1};
  132.    double A[] = { -0.334, 0.286 };
  133.    double X[] = { -0.14, -0.135 };
  134.    int incX = -1;
  135.    double Y[] = { 0.455, 0.358 };
  136.    int incY = -1;
  137.    double A_expected[] = { -0.264521, 0.0 };
  138.    cblas_zher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
  139.    {
  140.      int i;
  141.      for (i = 0; i < 1; i++) {
  142.        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher2(case 1455) real");
  143.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher2(case 1455) imag");
  144.      };
  145.    };
  146.   };
  147.  
  148.  
  149.   {
  150.    int order = 102;
  151.    int uplo = 121;
  152.    int N = 1;
  153.    int lda = 1;
  154.    double alpha[2] = {-0.3, 0.1};
  155.    double A[] = { -0.334, 0.286 };
  156.    double X[] = { -0.14, -0.135 };
  157.    int incX = -1;
  158.    double Y[] = { 0.455, 0.358 };
  159.    int incY = -1;
  160.    double A_expected[] = { -0.264521, 0.0 };
  161.    cblas_zher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
  162.    {
  163.      int i;
  164.      for (i = 0; i < 1; i++) {
  165.        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher2(case 1456) real");
  166.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher2(case 1456) imag");
  167.      };
  168.    };
  169.   };
  170.  
  171.  
  172.   {
  173.    int order = 102;
  174.    int uplo = 122;
  175.    int N = 1;
  176.    int lda = 1;
  177.    double alpha[2] = {-0.3, 0.1};
  178.    double A[] = { -0.334, 0.286 };
  179.    double X[] = { -0.14, -0.135 };
  180.    int incX = -1;
  181.    double Y[] = { 0.455, 0.358 };
  182.    int incY = -1;
  183.    double A_expected[] = { -0.264521, 0.0 };
  184.    cblas_zher2(order, uplo, N, alpha, X, incX, Y, incY, A, lda);
  185.    {
  186.      int i;
  187.      for (i = 0; i < 1; i++) {
  188.        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher2(case 1457) real");
  189.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher2(case 1457) imag");
  190.      };
  191.    };
  192.   };
  193.  
  194.  
  195. }
  196.